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o ■ 

^ ' Abstract. Infrared features of gluon propagator, ghost propagator, QCD running coupling and the 

Q , Kugo-Ojima parameter in lattice Landau gauge QCD are presented. The framework of PMS analysis 

■ suggests that there appear infrared, intermediate and ultraviolet regions specified by A^, j3o and 

j3i. The propagators and the running coupling of ^ > IGeV are fitted by the MOM scheme using 

■ the factorization scale as the scale parameter. The running coupling data of ^ < lAGeV are fitted by 
I the contour improved perturbation method. The Gribov copy problem is studied in SU{2), j3 = 2.2 
. samples by comparing the data gauge fixed by the parallel tempering method and the data gauge 

^ ' fixed by the straightforward gauge fixing. The Gribov noise effect turned out to be about 6%. 

^ : 
\o . 

5^ INTRODUCTION 

o : 

2 ! In Landau gauge QCD, it is necessary to restrict the gauge configuration expressed by the 

I link variable Ux,^ toQ.L = {U\ — dD{U) > , (9A = 0} which is called Gribov region (lo- 

■ cal minima) 1 2], where D is covariant derivative. Zwanziger argued that the uniqueness of 

^ \ the gauge field is guaranteed in its subset Al = {U\A=A{U),Fu{l) = Min^Ff/ {g) } , and 

g_) ' called the \l fundamental modular(FM) region[3]. Here the optimizing function Fu{g) 

1 

— ■ is defined in the case of t/-linear {A^^^ = -(t/x,^ - Ul^)\triess p.) and logC/ (t/^,^ = 

: ^-^) as, Fu(g) = L,,^ (1 -iRetrf/|^) , and Fu{g) = \\A8\\^ = L.,^tr (A^,,^Uf,^), 

■ respectively|4]. 

. • In Landau gauge, QCD running coupling can be extracted from three gluon coupling 

and/or ghost-gluon coupling. In terms of gluon dressing functiuon Z^(^^) and ghost 
dressing function G(^^)^which are measurable in lattice Landau gauge, renormalization 
group invariant quatity[5] (Xs{q^) = g^G{q^)^ZA{q^) /4-n can be measured. 

Colour confinement in infrared QCD is characterized by the Kugo and Ojima param- 
eter u{0) which is defined as 

x,y \ / xy P 

where u"^{p^) = 5"^u{p^), u{Qi) = —c. The sufficient condition of the colour confine- 
ment is m(0) = — 1 . We measure the running coupling and the Kugo-Ojima parameter on 
lattice and study the ambiguity due to Gribov copy. 



NUMERICAL RESULTS OF LATTICE LANDAU GAUGE QCD 



The gluon propagator and the ghost propagator 

The gluon propagator of colour SU {n) gauge is defined as 

= (1) 

and the gluon dressingf unction as Z^i^q^) = q^DA_{q^). Our result is shown in FigUJa) 
and is consistent with|0,[l]. 

The ghost propagator, which is the Fourier transform of an expectation value of the 
inverse Faddeev-Popov operator ^ = —dD = —d^{ \ —M) 

D"^{x,y) = (tr(A"x|(^[t/])-i|AM) (2) 

where the outmost () denotes average over samples U , is evaluated as follows. We 
take plane wave for the source bl^l = A e'^^ and get the solution of Poisson equation 
= (-A)-ibW. We calculate iteratively =M^['](x)(z = l,---,yt- 1). The 

iteration was continued until MajCx|^[^](x)|/MaJCx|Efj/ 0f'^(x)| < 0.001 ~ 0.01. The 
number of iteration k is of the order of 60, in SU{2), 16^ lattice, and of the order of 
100 in 5t/(3). We define <^*(x) = Z^=i ^^'^x) and evaluate (A«e'*^,4>*(x)) as the ghost 
propagator from colour b to colour a. We also used the conjugate gradient(CG) method 
and observed that the data agree except at the lowest momentum point of 48^^. Results 
of CG method are shown at this momentum point. 
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FIGURE 1. (a)The gluon dressing function as the function of the momentum q{GeV). )3 = 6.0, 
24'*(triangle), 32'^(diamond) and j3 =6.4, 48'*(star) in log[/ version. The fitted Hne is that of the MOM 
scheme. (b)The ghost propagator as the function of the momentum ^(GeV). )3 =6.0, 24"*, 32^* and /3 =6.4, 
48"* in logU version. The fitted hne is that of the MOM scheme. 

The gluon and the ghost dressing functions are fitted by using the framework of 
principle of minimal sensitivity(PMS)|'T^ or the effective charge method lfllll . in which 
the running coupling h{q) is parametrized as 



Kq) = W^){l+W^)'(i82//3o-(/3i//3o 



+ W^)'^(i83//3o-(/3i//3o)^) + ■■•}. (3) 



The parameter yj^{q) can be expressed by the parameter y defined as a solution of 



j8olog^ = - + ^log(/3o3;) (4) 



'A2 y j8o 

where A characterizes the scale of the system, and the function 



Kq\y) = - + ^iog{M - /3oiog(^VAy . (5) 
y Po 



In PMS, y is treated as a function of q^. However, in this work we fix the scale by the 
factorization scale /i = l.91GeV at which renormalizable quantity can be approximately 
factorizable into scheme independent and dependent parts lUll lsll. 

In order to be consistent with the MS scheme, we define the variable z as 

Z = _g(-l-fe/2c) ^ _}_(^_^yb/c^iK7t 

= -Z{q^)e^''^ (6) 



where t = \og{q^/k^), c = /3i//3o = 51/22, b = /3o/2 = 1 1/2, = {2c/b)-^-/'^A^, 

K = -b/lclM,^. When we fix y by the solution of eq.® with /i = 1.91GeV, (MOM 
scheme), we find that in the gluon dressing function the Landau pole at z = l/e remains, 
and another pole appears at z ~ 0. 17. When y is chosen as q^ dependent, the three regions 
< z < 0. 17, 0. 17 < z < l/e and l/e < z can be continuously connected! 12], but there 
is a subtle problem of PMS in low energv l.141 and leave the problem to bridge the three 
regions to the future. 

The gluon dressing function q^DA{q^) and the ghost propagator DQ{q^) in MOM 
scheme are plotted in FigCta) and (b), respectively. They are singular dX q = A^ ~ 
Q.lSGeV which should be washed away by the non-perturbative effects. 



The QCD running coupling and the Kugo-Ojima parameter 

The QCD running coupling as{q) turned out to have a peak of the order of 1 at 
q ~ O.SGeV and decreases to a finite value dX q = 0. We fitted the data by the contour 
improved perturbation series, which is a way to make a resummation of the series of 
coupling constant. 

The effective running coupling in the MS scheme is expressed by the series of cou- 
pling constant h^"'^ as[ 12]. 

= h^"\\ +Ai/zW +A2/z(")2 + ■ ■ • +A„/zW") (7) 



The result of MOM scheme using y = 0.01594 is shown by the solid line in Fig.|3a). 
The lattice data of 24"^, 32'^ and 48"^ and the MOM scheme agree in O.SGeV <q< 2GeV 
but slightly overestimates in ^ > 2GeV. 

In contour improved perturbation method, physical quantities M are expressed in a 
series physical quantities M are expressed in a series 



(- 



«=1 



-1 



)"dd 



(8) 



(9) 



" 27tJ-nc{\+W{Z{q^)e'^(^ 

Terms in the series (EJ) have alternating sign and A3 is not known. By choosing half of 
A2 in the series, we can obtain a good fit of the data(Fig|2tb)). 
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FIGURE 2. (a) The running coupling a,{q) of /3 = 6.0, 24'^(box), 324(triangle), j3 = 6.4, 324(diamond) 
and 48'*(star) as a function of momentum q{GeV) and the resuh of the MOM scheme, (b) The running 
coupling as{q) as a function of momentum q{GeV) of j3 = 6.4,48^ lattice. The solid line is the result of 
using — 231 MeV. Dotted hne is the result of e^^l^^A.jj^ and including half of A2. Dashed hne is 
the result of MOM scheme. 

The Kugo-Ojima parameter depends slightly on definitions, ?7— linear or \ogU and c 
becomes larger as the lattice size becomes large. However the value saturates at about 
0.8. 

TABLE 1. The Kugo-Ojima parameter c, trace e/d and horizon function 
deviation h in {/—linear and \ogU version. j3 — 6.0 and 6.4. 



/3 


L 


c\ 


exjd 


hi 


C2 




hi 


6.0 


16 


0.576(79) 


0.860(1) 


-0.28 


0.628(94) 


0.943(1) 


-0.32 


6.0 


24 


0.695(63) 


0.861(1) 


-0.17 


0.774(76) 


0.944(1) 


-0.17 


6.0 


32 


0.706(39) 


0.862(1) 


-0.15 


0.777(46) 


0.944(1) 


-0.16 


6.4 


32 


0.650(39) 


0.883(1) 


-0.23 


0.700(42) 


0.953(1) 


-0.25 


6.4 


48 








0.793(61) 


0.982(1) 


-0.19 



THE GRIBOV COPY PROBLEM 



We applied the parallel tempering(PT) method, which is used in finding global minimum 
in spin glass systems, to perform the FM gauge fixing of SU{2) /3 = 2.2, 16^ lattice 
configuration and compared with ensemble of 1st copy(Gribov copy obtained by the 
straightforward gauge fixing)! 6]. The Kugo-Ojima parameter c of PT becomes smaller 
than 1st copy by about 4%. Gribov noise in the ghost propagator is similar to Cucchieri's 
noise] 9] at /3 = 1.6, i.e. the infrared singularity of PT is weaker than 1st copy by about 
6%. The SU(2) gluon propagator, ghost propagator and running coupling in q > IGeV 
region are consistent with results of[15]. The running coupling in our simulation has a 
peak at around q = IGeV and is suppressed in the infrared. 

We aimed to detect in the lattice dynamics, the signal of Kugo-Ojima confinement cri- 
terion derived in the continuum theory, formulated in use of Faddeev-Popov Lagrangian 
and BRST symmetry. We also noted Zwanziger's horizon condition, based on the lattice 
formulation, coincides with Kugo-Ojima criterion||2l0]. However, our present data are 
not satisfactory to prove or disprove the confinement criterion. The colour off-diagonal 
antisymmetric part of ghost propagator lfl6[ \vh does not appear in Landau gauge, and 
the off-diagonal symmetric part has vanishing statistical average but has fluctuation pro- 
portional to {qa)^^. 
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